DESY 09-033 
SFB/CPP-09-26 
LPT-ORSAY 09-17 
HU-EP-09/12 

Lattice calculation of the Isgur-Wise 

g: functions T1/2 and r3/2 with dynamical quarks 

o ■ 

(N 

P5 ' Benoit Blossier'^'^, Marc Wagner'^, Olivier Pene^ 

in ! " DESY, Platanenallee 6, D-15738 Zeuthen, Germany 



^ Laboratoire de Physique Theorique (Bat. 210), Universite Paris-Sud XI, 
Centre d'Orsay, 91405 Orsay-Cedex, France. 



Oh; 

QJ \ Humboldt-Universitat zu Berlin, Institut fiir Physik, NewtonstraBe 15, D-12489 Berlin, 

Germany 



oo 
(N 

s 



^ . March 13, 2009 




Abstract 

We perform a dynamical lattice computation of the Isgur-Wise functions Ti/2 and at zero re- 
coil. We consider three different light quark masses corresponding to 300 MeV < mps < 450 MeV, 
which allow us to extrapolate our results to the physical u/d quark mass. We find 
''"1/2(1) = 0.296(26) and T3/2(1) = 0.526(23). Uraltsev's sum rule is saturated up to 80% by 
the ground state. We discuss implications regarding semileptonic decays B Xclv and the 
associated "1/2 versus 3/2" puzzle. 



1 Introduction 



The semileptonic decay of B mesons into positive parity charmed mesons (often referred to as 
D**'s) is an important and debated issue. Important, because no accurate measurement of the 
Vcb CKM angle will be possible, if these channels, which represent about one quarter of the 
semileptonic decays, are not well understood. Debated, because there seems to be a persistent 
discrepancy between claims from theory and from experiment [T]. 

Two types of Z)**'s are seen, two "narrow resonances" and a couple of "broad resonances", 
grossly speaking in the same mass region. While experiments point towards a dominance of 
the broad resonances in semileptonic decays, theory, when using the heavy quark limit, points 
rather towards a dominance of the narrow resonances. To clarify the situation ref. [1] called for 
actions on both the experimental and the theoretical side. 

The theoretical argument relies on a series of sum rules [21 [3] derived from QCD comforted 
by model calculations [H O [6]. Lattice calculations are needed to give a more quantitative 
prediction stemming directly from QCD. A preliminary computation was performed in [7], but 
only in quenched QCD and with a marginal signal-to-noise ratio. In this letter we report on 
the first unquenched computation using Nf = 2 flavor gauge configurations with Wilson twisted 
quarks generated by the European Twisted Mass Collaboration (ETMC). The spectrum of 
heavy-light mesons in the static limit has already been reported [101 Hi]- 

1.1 Spectrum in the heavy quark hmit 

We treat both b and c quarks via static Wilson lines, i.e. consider their infinite mass limit. In 
this limit the meson spectrum is constructed by combining the spin 1/2 of the heavy quark 
with the total angular momentum and parity of the light degrees of freedom (light quarks 
and gluons) [121 \T3\ I14j . The two lightest negative parity mesons B and B* (or D and D*) 
are degenerate and described by the same S = (l/2)~ state of light particles. The lightest 
(non-radially excited) positive parity states can be decomposed into two degenerate doublets: 
P_ = (1/2)+ and P+ = (3/2)^". The total angular momenta of the P_ (P+) mesons are O"*", 
1+ (1"*^, 2+). The mixing between the two l"*" states is suppressed in the heavy quark limit. 

It is generally believed that the narrow (broad) resonances are of the P^. (P-) type, since in the 
heavy quark limit they decay into D^*^tt via a D (S) wave. The D wave decays are supposed to 
be suppressed by a centrifugal barrier, if the final state momenta are not too large. 

1.2 Decay form factors in the heavy quark hmit 

In the heavy quark limit the semileptonic decay of a pseudoscalar meson into D** is governed 
by only two form factors [H], ti/2{w) and T3/2('w), where w = vb ■ vd** > 1 with vb and vd** 
denoting the four-velocity of heavy-light meson H being defined by vh = Ph/ttih- Uraltsev has 
proven the following sum rule [3]: 




(1) 



n 
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where t^^2^w) (t|"2(^))) ?^ = 0, . . . , oo are the form factors for the decay into the {P-) meson 
and the tower of its radial excitation^, w = 1 corresponds to the zero recoil situation, i.e. the B 
and the D** meson have the same velocity. Eqn. ([T|) is one of the major among many theoretical 
arguments in favor of the narrow resonance dominance [1]. 

Our goal in this paper is to make a direct lattice calculation of Ti/2(1) and t^/2{^) using static 
quarks represented by Wilson lines [15]. However, there is the problem that the B D** decay 
amplitude is suppressed at w = 1 due to vanishing kinematical factors, which multiply Tj(l). 
This is also a centrifugal barrier effect, i.e. it is impossible to give angular momentum to a meson 
at rest. Consequently, a computation of the weak current matrix element will trivially give zero. 
To overcome this difficulty, we use a method, which amounts to compute the operator matrix 
element based on an expression of the derivative of that matrix element in terms of the recoil 
four- velocity of the final meson [16\ [7]. Thanks to the translational invariance in time of the 
heavy quark Lagrangian this is then proportional to Tj{l){TnHj — mn), j = 1/2,3/2 (cf. eqns. 
(jlip and (|12p ). The mass splittings niH** — mn have already been computed in the static limit 
with precisely the same setup we are using in this paper [10^ 111] , i.e. by using Nf = 2 ETMC 
gauge configurations. We are thus in a position to compute ti/2(1) and t^/2{^) and to confront 
it with the Uraltsev and other sum rules as well as with other non-perturbative estimates (QCD 
sum rules, quark models). 

Our work should help to clarify the situation in the heavy quark limit. A fair comparison 
with experiment further needs to estimate the systematic error stemming from the heavy quark 
limit. After all, the charm quark is not so heavy. The authors of [5l[6] argue that large ©(l/mg) 
corrections are present. This issue can also be addressed by lattice QCD, but in this work we 
restrict our computations to the static limit. 

The paper is organized as follows. In section [2| we recall the method used to compute ti/2(1) 
and T3/2(l)- In section [3] we report on the lattice calculation of ri/2(l) and T3/2(1)- In section 
m we perturbatively compute the renormalization constant of the heavy-heavy current and we 
conclude in section [5j 



2 Principle of the calculation 

To compute the zero-recoil Isgur-Wise functions ri/2(l) and T3/2(1) by means of lattice QCD, 
we use a method proposed in [7]. We remind it here just for comfort of the reader. 

The method consists in using a series of relations derived in ref. [16]. With v' = (1,0,0,0) and 
V = v' + v±_ denoting the velocities of the ingoing and outgoing mesons, where v± is spatial up 
to higher orders in the difference v' — v, we assume that for some Dirac matrix F/ 

{H**{v')\Q{v')TiQ{v)\H^*Hv)) = tTv^mTjiw) + ■ ■ ■ (2) 

Here w = v ■ v', j = 1/2, 3/2 and l,m = 1,2,3 are spatial indices, tf^ is a tensor, which 
depends on the final state (H**) and the initial state {H* or H), and Q{v) is the static quark 
field in Heavy Quark Effective Theory. The dots represent higher order terms in v' — v. From 

^By definition Tj(w) = t^°\w), j = 1/2, 3/2. 
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translational invariance in time direction 

-ido{H**{v')\Q{v')TiQ{v)\H(*^{v)) 

= -i{H**{v')\Q{v') \Ti + v 'l 

= tY'v±rnTj{w)(^mH" - TTlH^ H 



Q{v)\H'~*\v)) 



Then we use the field equation {v ■ D)Q{v) = 0: 
D'^Qiv') = , D^Q{v) = -{D-v±)Q{v), 
whence from eqn. ^ 

i{H**{v')\Q{v')Ti{D ■ v^)Q{v)\H'^*Hv)) = tr^'WiH (m^- - mj^) + • • • 
which, in the limit v± 0, converges to the relation 



i{H**{v)\Q{v)TiD"^Q{v)\H^*Hv)) = trTj{l)(mH** - mn 



Applying eqn. ([2]) to the J = Hq state we get from ref. [T7j 

{H^{v')\A\H{v)) = -T^,^{w)vu, 



(3) 



(4) 



(5) 



(6) 



(7) 



where Ai is the axial current in spatial direction i, and where the normalization of the states is 
\/\/2m times the one used in ref. [T7]. From eqn. ([7]) follows 



{HQ{v)\AiDj\H{v)) = igij(mH'^-mH)Ti/2{l). 
Analogously for the J = 2 H2 state we have 
{H*2{v')\A,\H{v)) ^ V3r3/2M6f^;±, + ---, 
where £*■' is the polarization tensor, whence 
{Hl{v)\AiD,\H{v)) = -iV3(mH^ - MH)Ty2{l)e*y 

Finally t^/2(1) ^^'^ '^3/2(1) can be obtained from the following matrix elements: 
{H*\Qj5lzD,Q\H) 



(8) 



(9) 



(10) 



n/2(l) 
7-3/2(1) 



rriH* - ruH 



{m\Qj5{lxD,-lvDy)Q\H) 



(11) 
(12) 
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There is no mixing of the operators AiDj with dimension 3 (hence hnearly divergent) heavy- 
heavy operators to be feared on the lattice: indeed we are interested in a parity-changing 
transition and ah dimension 3 operators have vanishing matrix elements between positive and 
negative parity state^. There are no logarithmic divergence either thanks to the vanishing of the 
vector and axial currents' anomalous dimension in HQET at zero recoil. By consequence there 
is no conceptual issue concerning the extrapolation to the continuum limit of such a calculation. 
It needs only a finite renormalization constant to match the lattice result with a continuum-like 
scheme value, as we will discuss in Section HI 



3 Lattice computation of ri/2 and at zero recoil 
3.1 Simulation setup 

We use Nf = 2 flavor 24^ x 48 Wilson twisted mass gauge configurations produced by the 
European Twisted Mass Collaboration (ETMC). Here we only give a brief summary of the 
setup, which is explained in detail in [T9 l [20]. 

The gauge action is tree-level Symanzik improved [21] with /? = 3.9 corresponding to a lattice 
spacing a = 0.0855(5) fm: 



Sg[U] = 




where = 1 — 8bi and 6i = —1/12. 

The fermionic action is Wilson twisted mass with two degenerate flavors |22^ [23] : 
Sf[x,X,U] = a''^x(a;)(^W + ^/^q75T3)x(3;), (14) 

X 

where 

= ^(7m(Vm + V;) -aV*V^) +mo, (15) 

and V* are the standard gauge covariant forward and backward derivatives, uiq and are 
the bare untwisted and twisted quark masses and x = (x^"^ > X^^^) a-rs the fermionic fields in the 
twisted basis. 

We consider three different values of the light quark mass, which amount to "pion masses" in 
the range 300 MeV < mps ^ 450 MeV (cf. Table [1]). mo has been tuned to its critical value at the 
lightest /Xq value, i.e. at fiq = 0.0040. 

^Of course the situation is different by instance for the matrix element {H\hY)'^h\H), related to the HQET 
parameter Ai or the kinetic momentum for which a subtraction is necessary to its computation on the lattice 
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mps in MeV 


number of gauge configurations 


0.0040 


314(2) 


1400 


0.0064 


391(1) 


1450 


0.0085 


448(1) 


1350 



Table 1: twisted quark masses fiq, pion masses mps and number of gauge configurations. 



3.2 Static and light quark propagators 

The propagator of a static quark is essentially a Wilson line in time direction: 
[Qix)Qiy))^_ = 5(^\^-y)U('^^^^\x;y)(Q{yo-xo)^^ + Q{xo-yo)^^y (16) 

where (. . ■)q^q denotes the integration over the static quark field and U{x;y) is a path ordered 
product of links along the straight path from x to y. To improve the signal-to-noise ratio we 
use the HYP2 static action [25l [26| [27] . 

For the light quarks we use four stochastic spin diluted timeslice propagators {Z2 x Z2 sources 
with randomly chosen components ±1 it i) for each gauge configuration. For details we refer to 
[TT] . where exactly the same setup has been used. 

3.3 Static-light meson creation operators 

In the static limit there are no interactions involving the heavy quark spin. Therefore, it is 
convenient to classify static-light mesons according to , where j denotes the angular momen- 
tum of the light degrees of freedom and V parity. In particular we are interested in the sectors 
= (l/2)~, = (1/2)"'" and = (3/2)"*". We label the corresponding static-light mesons, 
i.e. the ground states in these angular momentum/parity sectors, by S, P- and P+ respectively. 

To create such static-light mesons on the lattice we use operators 

C)(r)(x) = Q(x) r(n)[/(x;x + rn)x^''^(x + rn), (17) 

n=±ei,±e2,±e3 

where Q creates a static antiquark at position x, x^"^ creates a light quark in the twisted basis 
at position x -|- rn, [/ is a product of spatial links along the straight path between x and x -|- rn, 
and r is a combination of spherical harmonics and 7 matrices yielding a well defined behavior 
under cubic rotations (cf. Table [2]). 

To optimize the ground state overlap of these static-light meson states, we use Gaussian smearing 
|28j for light quark operators and APE smearing ^29j for spatial links (parameters Kcauss = 0.5, 
A^Gauss = 30, OAPE = 0.5, iVAPE = 10 and r = 3 as in [IT])- 
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r(n) 


Oh 


3 


75 


Ai 


1/2 , 


7/2 , ... 


1 




1/2 , 


7/2 , ... 


Ixf^x - lyfiy (and cyclic) 


E 


3/2 , 


5/2 , ... 


Ibilxfix - lyTiy) (and cyclic) 




3/2 , 


5/2 , ... 



Table 2: static-light meson creation operators. 



3.4 Static-light meson masses 

Since we work in the twisted basis, where each of the operators listed in Table [2] creates both 
V = + and V = — states, it is convenient to determine V = + and V = — static-light meson 
masses from the same correlation matrix. 

For 5 and P- we compute the 2x2 matrix 

CjK{t) = U0^^-'\t)]^0^^^\Q)), (18) 



where Tj G {75 , 1}, and solve the generalized eigenvalue problem 

CjK{t)v^;^\t) = CjK{to)vP{t)X^''Ht,to) , to = 1 (19) 

(cf. [30^ I31j). The meson masses m{S) and m{P^) are determined by performing minimizing 
fits to effective mass plateaus, 

it) - Inf^^^^ML^ (20) 

at large temporal separations t (as indicated in Figure [1] our fitting range is 6 < t < 11). 
The parity of the corresponding states, i.e. whether it is S or P_, can be extracted from the 
eigenvectors Uj"^ (for a detailed discussion, of how to identify parity, cf. [H]). Results of meson 
masses and mass differences and corresponding reduced values are listed in Table [3l 

For m{Pj^) we proceed analogously this time computing the 2x2 matrix (jl9p . where 

rj G {^xnx - lyfly , Ibilx^x " 7y^?/)}- 

By solving the generalized eigenvalue problem (|19|) we have also obtained appropriate linear 
combinations of twisted basis meson creation operators with well defined parity. To be more 
precise the operators 

O(^) = t;(,f (t)0(T5) +r;f)(t)0« (21) 
O'^P-) = i;(f-)(t)0(^^)+t;f-)(t)0« (22) 

(n{P+) = 7/"^+^ ('+^C)(7a:na:-7^;'^!;) I (f\(nil5i-yxn:,-'yyny)) /2o^ 
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(iq = 0.0040 - ma„aai„a(S) and mo„j,ai„a(PJ 



|iq = 0.0040 - mg„^„,i„^(PJ 
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(iq = 0.0064 - ma„aai„a(S) and mo„j,ai„a(PJ 
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Figure 1: effective masses for S, P_ and P+ for /iq G {0.0040, 0.0064, 0.0085}. 



create static-light meson states, which have the same quantum numbers j as the states of 
interest, \S), and \P+) respectively. Since the t dependence of the eigenvectors Uj"^ is very 
weak [32] , results are essentially unaffected by the choice of t (we have used t = 6 for all results 
presented in the following). 
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m{S) 


xVdof 


m(P_) 


xVdof 




xVdof 


0.0040 


0.3987(19) 


1.79 


0.5670(60) 


1.69 


0.6101(66) 


2.46 


0.0064 


0.4061(17) 


1.93 


0.5877(67) 


0.45 


0.6121(64) 


3.01 


0.0085 


0.4104(17) 


2.23 


0.6095(65) 


0.49 


0.6283(41) 


0.87 





m{P-) — m{S) 


m{P+) — m{S) 


0.0040 
0.0064 
0.0085 


0.1683(65) 
0.1817(69) 
0.1991(63) 


0.2114(62) 
0.2060(63) 
0.2179(41) 



Table 3: static-light meson masses and mass differences for /Xq G {0.0040, 0.0064, 0.0085}. 



3.5 Two-point functions and their ground state norms 

After having obtained the linear combinations of twisted basis operators ()2ip to ()23p the two- 
point functions 

((O(^)(t))^0(^)(0)) , ((o(^-)(t))^O(^-)(0)) , ((o(^+)(i))^O(^+)(0)) (24) 
are straightforward to compute. 

From these two-point functions we also determine the ground state norms of the corresponding 
sectors, N{S), N{P^) and N{P^), by fitting exponentials at large temporal separations. To 
be more precise, we obtain e.g. N{S) by fitting iV(5)2e-'"* to {(O^^^t))^ O^^\0)) with N{S) 
and m as degrees of freedom. Results and corresponding reduced values are listed in Table H] 
(fitting range 6 < t < 12). 



A'q 


N{S) 


xVdof 


iV(P_) 


xVdof 


N{P+) 


xVdof 


0.0040 


0.3271(26) 


0.21 


0.2998(93) 


0.33 


0.1139(26) 


1.43 


0.0064 


0.3358(20) 


0.23 


0.3074(87) 


0.13 


0.1120(27) 


1.68 


0.0085 


0.3397(22) 


0.22 


0.3139(103) 


0.08 


0.1212(22) 


0.28 



Table 4: ground state norms for ^Uq G {0.0040, 0.0064, 0.0085}. 



3.6 Three-point functions and form factors ti/2 and 

In analogy to effective masses we define effective form factors 

''"l/2,efrcctivc(*0 " ^2) 

1 NiP^) N{S) {[o(P-\to)y {Qj,^,D,Q)it,) 0(3\t2)) 

(m{P.)-m{S)) ((o(^-)(io))^0(^-)(ti)) {{oi^Kh))^ 0(S) (t2)) 
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''"3/2,cfrcctivc(*0 - h,ti - t2) 



1 



N{P+) N{S) ((o(^+)(to)) {Ql5{lxD,-7yDy)Q){h) 0(^\t2) 



V6 (m(P+)-m(5)) ((o(^+)(io))W)(ti)) ( (o(^)(ti))^0(S)(t2 



(26) 



{Zx> = 0.976 is a lattice renormalization constant, which we derive and discuss in detail in 
section m. These effective form factors are related to ti/2 and r3/2 via (fTT|) and p^ : 



'^1/2(1) — , , , '^l/2,cffcctivc(*0 - - ^2) 

'^3/2(1) = ^ ^ . 7-3/2.efrectivc(*0 - - ^2)- 



(27) 
(28) 



Computation of the three-point functions appearing in (j25p and (j26p is again straightforward. 
We chose to represent the covariant derivative acting on the static quark field symmetrically by 



DjQ{^,t) 



-([/j(x,t)Q(x + e„t) 



t 



Uj (x -ej,t)] (5(x - Bj , t) 



(29) 



To optimally exploit our gauge configurations and propagator inversions, we average over all 
three-point functions, which are related by the lattice symmetries 75 hermiticity, parity, time 
reversal, charge conjugation and cubic rotations. 

The resulting effective form factors Tx/2,cffcctivc(^o ~ ^ii *i ~ ^2) and T3/2,cfrcctivo(*o — — ^2) are 
shown in Figure [2] as functions of to — ti for fixed to — t2 £ {10, 12}. Within statistical errors 
these effective form factors exhibit plateaus for to — ti ^ {to — i2)/2, i.e. when both temporal 
separations, to—ti and ti—t2, are large. We determine T1/2 and T^/2 by performing minimizing 
fits to the central three data points as indicated in Figure [2j Results for to ~ ^2 = 10 and for 
to — 12 = 12, which are listed in Table EJ are in agreement within statistical errors. We consider 
this a strong indication that contributions from excited states at these temporal separations are 
essentially negligible and that the plateaus of the effective form factors indeed correspond to T1/2 
and T3/2- In the following discussions we only quote the numbers corresponding to to — t2 = 10, 
since their statistical errors are significantly smaller than those for to — t2 = 12. 





to — t2 


n/2 


■^3/2 


T3/2/T1/2 


(.r3/2? - {ti/2? 


0.0040 


10 


0.299(14) 


0.519(13) 


1.74(9) 


0.180(16) 




12 


0.267(26) 


0.536(25) 


2.01(21) 


0.216(30) 


0.0064 


10 


0.312(10) 


0.538(13) 


1.73(6) 


0.193(13) 




12 


0.278(19) 


0.549(21) 


1.98(14) 


0.225(23) 


0.0085 


10 


0.308(12) 


0.522(8) 


1.69(6) 


0.177(9) 




12 


0.287(24) 


0.544(14) 


1.90(17) 


0.214(21) 



Table 5: ri/2 and r3/2 for to - ^2 G {10 , 12} and G {0.0040 , 0.0064 , 0.0085}. 
As expected from operator product expansion, r3/2(l) is significantly larger than ri/2(l)- More- 



over the Uraltsev sum rule [3], 



1 

4' 



(30) 



is almost fulfilled by the ground state contributions Ti/2i^) = ''"1/2(1) and ''"3/2(1) = T3/2(1)- 

Finally we use our results at three different light quark masses (cf. Table [T|) to perform a linear 
extrapolation of the form factors in (mps)^ to the physical u/d quark mass (mps = 135 MeV). 
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Figure 2: effective form factors Ti/2,cfrcctive and r3/2,efrective ^ov to — *2 G {10 , 12} and 
/iq G {0.0040, 0.0064, 0.0085}. 
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Results are shown in Figure [3] and Table [6l The qualitative picture for u /d quark masses is 

'3/2 



the same as for the heavier masses used directly in our simulations: tT/^^^^''{1) = 0.526(23) is 



significantly larger than t^^^^''{1) = 0.296(26) supporting the "theory expectation" that a decay 
of a -B meson to a j = 3/2 P wave D meson is more likely than to a j = 1/2 P wave D meson. 



extrapolation to ttie ptiysical u/d quark mass (tg - 12 = 1 0) 



0.6 

0.5 
0.4 - 
0.3 
0.2 
0.1 
- 



lattice results — ^— 
extrapolations =hii#ii 



0.1 0.15 
(mps)^ in GeV^ 



extrapolation to the physical u/d quark mass (Iq - f2 = 12) 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




lattice results — 
extrapolations aiii^iifc 



0.1 0.15 
(mpg)^ in GeV^ 



Figure 3: linear extrapolation of T1/2 and T3/2 to the u/d quark mass for to — *2 £ {10 , 12}. 



to — t2 


n/2(l) 


xVdof 


^■3/2(1) 


xVdof 


10 


0.296(26) 


0.34 


0.526(23) 


1.43 


12 


0.251(48) 


0.00 


0.536(43) 


0.12 



Table 6: linear extrapolation of T1/2 and t^/2 to the u/d quark mass for to — t2 € {10 , 12}. 



4 Perturbative renormalization of the static current Qj^jiDjQ 

In this section we derive the analytical formulae and give the numerical values of the renor- 
malization constant of the dimension 4 current Oij = Q^5^iDjQ computed at first order 
of perturbation theory for the HYP smeared static quark action and both the standard Wilson 
plaquette and the tree-level Symanzik improved gauge action. 



4.1 Definitions 

The bare propagator of a static quark on the lattice is 

a 



I _ g-ip4a _^ a^jji aT,{p) 
^ Z2hS''ip). 
Choosing the renormalization conditions 



ip4—fO 



ip4 , 6m 



-^{PA = 0) 
11 



a{6m + S(p)) 



(31) 



(32) 



implies 

The bare vertex function Vj^{p) is defined as 

V^ip) = {S''r\p)Y,e'^^^-y^{Q''{x)Of^m^{y)){S^)-\p) 

= ^{S'')-\p)Y.e^^^^-y^{Q^{x)Of^mHy)){S'')-\p\ (34) 

x,y 

where 

0,^(0) = ZvOf^iQ). (35) 
V^{p) can be written as 

^if(p) = {l + 5V)u{p)-ia^pju{p) ^ (l + 5V)y,f(p). (36) 
bV is given by all the IPX one-loop diagrams containing the vertex. 

4.2 Analytical formulae and results 

The notations used in this section and the Feynman rules are listed in appendix \F\ They are 
the same as in [33] except for the gluon propagator having the form 

D^, = C^'Df^-^ + A^, (37) 

[31], where Cq = cq + 8ci + 16c2 + 8C3 = 1, ci = — 1/12, C2 = C3 = for the case of the tree- level 
Symanzik improved gauge action and 

A^^ = 6f,^K^ + AL^^N^N^. (38) 

Finally Kf^ and L^i^ are complicated expressions, which do not need to be reproduced here. 
The only relevant features for this work are that A^i, is regular in the infrared regime and 

The static quark self-energy expressed at the first order of perturbation theory is given by 
S(p) = — (-Fi -|- F2), where Fi and F2 correspond to the diagrams shown in Figure Hl^a) and (b): 

A C -i(fc4+2ap4) 

Fi = 9o / h^hiiDii rjr— — ^ = Ff^'^ + Fi (39) 



12 



^plaq 



3a 



9o 



, 2W + a^A^ 1 - e-«(fc4+ap4) + g 



4 2 /■i^lH£;) + ELiG'l,H^) 1 



4EVl + 

^3 ^2 



,42 f z?i(-ij^) + Etigi»(-^-g) 



1-e^' 
1 



(40) 



3a' 



-i(fc4+2ap4) 



<7o / ^4i/l4,Ai,-^— — -^^^^^^^ 



-ap4^0 



M4 - i A^4 



3a ^ 7fc 2iiV4 + eM4 



+ 50*^4 
-3 -Ik 



Mi-iNi _ I ( Mi- iN4 

£ 

1 



2iNi + eM4 2 ^2iiV4 + eM4 



MlK + 'i{plK^ + nIK 



(41) 



iV|A = 4 ( D|A^|(^i + £44) + \Y, GliiK"" + miK';) + 2D4iV4 ^ G4^^^^L, 

^ i i=l 

3 \ 

+2 ^ G4iG4jiV,Ar.L,^.j 



(42) 



F2 = -l^^"'""' I h4^h4jD,, 



(43) 



pplaq 
^2 



2 3a 



2W 



-ap4^0 



2 3 



(l/a - .^4) j[ 



3 /^2 

4i 



2W^ 



(44) 



F' 



(45) 




P p+k P 
(a): sunset diagram 



P P 
(b): tadpole diagram 



Figure 4: self-energy corrections. 
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The factor 1/2 has been introduced to compensate the over-counting of the factor 2 in the 
Feynman rule of the two-gluon vertex, when a closed gluonic loop is computed. 

The other terms entering the above integrals cancel, because the contour can be closed in the 
complex plane without including the pole = —p^ + i ln(l + e). Finally we can write 



Fi = 



127r2 



(46) 



Fo = 



127r2 



l/a-ip4)(/|'">0 + /3(«.,Q))- 



(47) 



The linearly divergent part in 1 /a of the self-energy is given by 



So(aj) 



9l 



-o-o(ai) , o-Q = fi + fi + h + fs, 



127r2a 

while the wave function renormalization Z2h reads 



(48) 



Z2h{o:i) 



1 + 



127r2 



Z2 



/3 + /3-(/2 + /2) 



(49) 



In Table [7] we have collected the numerical values of /j, f-, ctq and Z2 for different kinds of static 
quark and gluonic actions. 

The vertex function V.f is obtained by writing 



B 



yO _L yl I y2 



(50) 



corresponding to the diagrams (a), (b) and (c) in Figure [5l The contribution is given by 
computing 



3a^ 



9o / hikhAiDkis\ii{k + ap)j 



g-i(fc4+2ap4) 



(1 _ g-j(A:4+ap4) _|_ g)2 



T/0,plaq I 

^3 ^ ^3 



(51) 



0,plaq 



ii 2 [ Dl + E ii 

3a ° Jk 2W + a2A2 



-i(fc4+2ap4) 



(]_ _ g-«(fc4+ap4) _|_ 




(a) 





Figure 5: operator corrections. 
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1 

X 



3'^°^^ A 2VF + a2A2 (2iiV4 + eM4)2 ^^^^ 



= ~^9o hikhiiAki (Vj + apj cos{kj)^ (1 - iap4) 

= -J91P, I {dIk^ + n!a) cos(fe,) ^^^^^ j = l,oV,/^Acos(fc,) 
The "sail diagram" has the following expression: 

= A^2 / cos ( ^ + apj ^ ^ , = T/^'Pl"^ + V,'i 



3a''J,2W + a^X^\ ' 2 7 2i sin (^Ei) + e^^e 

_£ 2 [ G^jNj 1 _ 2 2. /- g4,^i 

3a^°^^ ifc 2H^ + a2A2 2iiV4 + eM4 3^°'^^ A 2W^ + a2A2 ' 

G4j = A^4Glj- (53) 



2f sin (tePi 



= ^,gj[/.4.A,,(M,-ap,iV,)(l-z^) 

= '^glipj I N][AD4Uj + M^ , N^NjM^ = Y^^a^^r (54) 
•^'^ 1=1 

Note that the contribution of the sail diagram to the final result must be doubled, because 
the gluon leg can be attached to the static line in two different ways. Eventually the tadpole 
diagram is given by 

^/(«*) = -^3^5oPiy^^44 = -^Pj(/3(«i = 0) + /^(a. = 0,q)). (55) 
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We finally have 



{H**\Ofj\H) 



where 



-{H**\Of^\H){ai), 



Zv{ai 



Z2h{ai)[l + 5V{a{)^ 

2 



I.e. 



Zv{ai) 



Z2 + f4 + fi- 



(56) 

(57) 
(58) 

(59) 



The numerical values of are collected in Table [7] for the different kinds of static quark and 
gluonic actions. With the bare coupling = 6//3, the tree- level Symanzik improved gauge 
action at /3 = 3.9 and the HYP2 static quark action used in our simulations we obtain 
Zi5(tlSym,HYP2) = 0.976. 





ai = 


HYPl 


HYP2 


fl 


7.72 


1.64 


-1.76 


/((tlSym) 


2.10 


0.14 


0.83 




-12.25 


1.60 


9.58 


/^(tlSym) 


-3.43 


-0.12 


-1.50 


h 


12.23 


4.12 


5.96 


/^(tlSym) 


-2.10 


-0.14 


-0.83 


h 


-12.68 


-4.95 


-0.56 


/l(tlSym) 


3.56 


2.04 


1.67 




19.95 


5.76 


4.20 


Z2(plaq) 


24.48 


2.52 


-3.62 


Z2(tlSym) 


25.81 


2.50 


-2.96 


Zd{p\aq) 


11.80 


-2.43 


-4.19 


Zrf(tlSym) 


16.69 


-0.41 


-1.85 



Table 7: numerical values of the constants /i, /(, /2, /2, /s, f^, fi, f!i, cq, Z2 and z^ defined in 
the text; = denotes the unsmeared Eichten-Hill static quark action, while HYPl and HYP2 
are defined in [25] and [27] respectively; "plaq" denotes the standard Wilson plaquette gauge 
action, while "tlSym" denotes the tree-level Symanzik improved gauge action. 



5 Conclusions 

We have computed the form factors ti/2(1) and T3/2(1) in the static limit, which describe (in 
this limit) the decay B — > D**. This decay is presently a puzzle in the sense that sum rules 
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derived from QCD point towards a dominance of T3/2(1), while experimental indications point 
rather in the opposite direction. The aim of this paper has been to check the dominance of 
"^3/2(1) ™ ^ quantitative way. 

Our final result extrapolated to the physical u/d quark mass is given in Table [H Since we see 
no systematic dependence on the temporal separation — 12 except for an increase in statistical 
uncertainty, we keep the result at to — t2 = 10. To the statistical error we add a systematical 
error of 3% to account for the uncertainty in the computation of the renormalization constant 
Zx>, which was computed perturbatively. We make the "guesstimate" of 100% uncertainty on 
1 — Zj), which turns out to be very small. Notice that this uncertainty does not apply to the 
ratio T3/2(l)/''"i/2(l) both having the same Zx> (cf. eqns. (p5]) and (p6]) ). We have at this stage 
no way to estimate systematic uncertainties arising from finite lattice spacing and from finite 
volume. Therefore, we must consider the errors we quote as incomplete. We end up with 

ri/2(l) = 0.296(26) , rg/aCl) = 0.526(23) (60) 
^•3/2(1) 

Tl/2(1) 



1.6. ..l.J 



2 



2 



^3/2(1) - T3/2(l) ~ 0.17... 0.21 (61) 



in fair agreement with the qualitative claim that is significantly larger than Ti/2- Note also 
that Uraltsev's sum rule is almost saturated by the ground state contributions providing ~ 80% 
of the required 1/4 (cf. eqn. ([1])). 

This result does not differ qualitatively from the preliminary quenched computation [7]: 
ri/2 = 0.38(5) and T3/2 = 0.53(8). However, we consider the result presented in this paper as 
standing on a much firmer ground, because it is unquenched, and because the signal is much 
clearer and more stable thanks to better analysis procedures. Our result ()60p is also similar to 
the prediction of a Bakamjian-Thomas relativistic quark model P], when using a Godfrey- Isgur 
interquark potential: T1/2 = 0.22 and T3/2 = 0.54. 

Assuming that the heavy quark limit provides reliable indications and that the standard iden- 
tification of narrow D** resonances is correct (i.e. -Di(2420) (J = 1) and 1)2(2460) (J = 2) 
correspond to j = 3/2 mesons) this points towards the expected dominance of the semileptonic 
decay of B mesons into these j = 3/2 states over the decay into j = 1/2 states. The latter, 
labeled as Dq ( J = 0) and D[ ( J = 1) are identified to some broad structures, which are seen in 
the semileptonic B decay around similar masses (2200 MeV to 2600 MeV). Remember, however, 
that the predicted ratio of branching fractions Br(i? — > L'3*2)/Br(i? ^1*2) niainly governed 
by (t3/2(1)/ti/2(1))^ times a rather large ratio of phase-space factors. 

It is usually claimed from experiment that the decay into these broad resonances are not sub- 
dominant as compared to the narrow resonances. A recent analysis by BABAR \35\ [36] finds 
significant B — > D'^*-' it I u, but does not give the relative yield of narrow and broad resonances. 
In a recent paper by BELLE [37] the four D** states are distinguished. The B — > Dq I v is 
observed with a comparatively large signal and, assuming the heavy quark limit to be applica- 
ble, they fit T3/2(l) = 0.75 and Ti/2(1) = 1-28. Compared to our result ()60p this calls for two 
comments: 

(1) The T3/2(l) shows fair agreement between theory and experiment. This is encouraging, 
since the narrow resonances are experimentally rather well under control, i.e. the narrow 
resonances are well seen. 
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(2) The experimental Ti/2(1) is much larger than our prediction. Note, however, that BELLE 
does not see the other member of the j = 1/2 doublet, B — > D'^lv. This is puzzling and 
the discrepancy concerning T2/2(1) should not be taken as final. 

In view of the impressive convergence of almost all theoretical estimates o/t]^/2(1) '^'^^ "^3/2(1); 
in view of our confidence that the result presented in this paper stands on a firm ground, we 
believe that one can consider as established that QCD predicts a clear dominance of the decay 
into j = 3/2 in the static limit. 

It still remains to be solved, how to saturate the inclusive semileptonic branching ratio, in 
other words what to add to the B D^*^ I v and to the narrow D** resonances. The analyses 
performed on Class I non-leptonic B — > D**7r decay do not find any trace of broad structures 
[381 I39| . Invoking factorization, theoretically well under control for this kind of process this 
naturally leads again to ti/2(1) < "^3/2(1)- 

Experimental work still has to be done. On the theory side, beyond doing the computation 
at another finer lattice spacing to be able to perform a continuum extrapolation (theoretically 
well defined, as recalled in Section [2]), an estimate of the l/rn-c corrections would help a lot. To 
explore that issue a promising method used to study the B — > / v form factors at non-zero 
recoil [401 [41] might be helpful. The contributions of other states such as negative parity radial 
excitations should also be considered. 

Let us conclude by insisting that the issue at clue is of important relevance: any accurate 
estimate of the Vch parameter of the standard model will never be fully convincing as long as 
the "1/2 versus 3/2 puzzle" remains unsolved. 
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A Feynman rules 

The lattice HQET action is 

5HQET = ^ (Q\n) (Q{n) - uf^^^'in - 4)Q(n - 4)) + abmQ\n)Q{n)\ , (62) 
where UY^^ in) is a link built from hypercubic blocking. 
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We will use in the rest of this appendix the following notations taken from \42\ 

Jp J^n/a (27r)* Jp- J_^/a (27r)'^ ^ 





-1; 






Tx = 


sin(aA;_ 




/a 


Cm = 


cos 1 — 


= 


COS 1 ■ 



= ^Ar2 + ^ , ^' = 2argsh(^). 



In Fourier space the action at 0(50) is given by 

^HQET ^ /■igt(p)(i_e-iP4a)Q(p)+^,„gt(p)g(p) 
Jp a 

+i90 [ [ [5{q+p'-p)Q^p)Bt{q)T''Q{p')e-'^P'+P'^^"^ 

J p J p' J q 

+ ^ / IJ l^'l + r+p' -p)QHp)Bt{q)Bl{r)T-T'Q{p')e-'^P^+P''^^^^ 



ip Jp' J q J r 

The block gauge fields B'^ can be expressed in terms of the usual gauge fields: 



00 

n) 



n=l 



where -B^"^ contains n factors of A. At next to leading order, it was shown that we only 
^(1) 



(A:) = Y.^Ak)Mk) 
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6^,uD^{k) + (1 - V)G'm.(^) 



Dp{k) 



di = -01(1 + 02(1 + 03; 



N.Njdi-d2 



Ni + N, 



^+d: 



-0102(1 + 203 



801O2O3. 



(77) 
(78) 

(79) 
(80) 



The Feynman rules are the following: 



heavy quark propagator 



vertex V, 



p,hhg 



vertex V/^l;f^f^gg{p,p') 
gluon propagator in the Feynman gauge 



a(l 



~tp4 



\agl5pA5yA Y.p,u hppKaiT", T^}e 
a'(Co-^ V5"'(2^ + a'A2)-i + A^,) 



Note that p' and p are the in-going and the out-going fermion momenta, respectively. We also 
introduce an infrared regulator A for the gluon propagator. We symmetrize the vertex V^^^^j^gg by 
introducing the anti-commutator of the SU{2>) generators normalized by a factor 1/2. The gluon 
propagator and the vertices are defined with the A field. At one-loop the infrared regulator to the 
gluon propagator that we have chosen is legitimate, because no three-gluon vertex is involved. 
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